R = MDOMReflection(Layer(Material('Al',2000),Inf)); R.N_in = 50; R.N = 151; R.Calculate;
for i=1:size(R.Rm,3)
Cn(i,1) = interp1(R.x,R.Rm(1,:,i),cosd(25));
end
R = MDOMReflection(Layer(Material('Al',1180),Inf)); R.N_in = 50; R.N = 151; R.Calculate;
for i=1:size(R.Rm,3)
Cn(i,2) = interp1(R.x,R.Rm(1,:,i),cosd(25));
end
R = MDOMReflection(Layer(Material('Al',500),Inf)); R.N_in = 50; R.N = 151; R.Calculate;
for i=1:size(R.Rm,3)
Cn(i,3) = interp1(R.x,R.Rm(1,:,i),cosd(25));
end
E0 = [2000 1180 500];

plot(Cn);